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In the proximity of melting transitions of artificial and biological membranes fluctuations in enthalpy, area, volume and 
concentration are enhanced. This results in domain formation, changes of the elastic constants, changes in permeability and 
slowing down of relaxation processes. In this study we used pressure perturbation calorimetry to investigate the relaxation 
time scale after a jump into the melting transition regime of artificial lipid membranes. This time corresponds to the 
characteristic rate of domain growth. The studies were performed on single-component large unilamellar and multilamellar 
vesicle systems with and without the addition of small molecules such as general anesthetics, neurotransmitters and 
antibiotics. These drugs interact with membranes and affect melting points and profiles. In all systems we found that heat 
capacity and relaxation times are related to each other in a simple manner. The maximum relaxation time depends on the 
cooperativity of the heat capacity profile and decreases with a broadening of the transition. For this reason the influence of a 
drug on the time scale of domain formation processes can be understood on the basis of their influence on the heat capacity 
profile. This allows estimations of the time scale of domain formation processes in biological membranes. 
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Introduction ations in the melting regime of membranes are accompanied 

by changes in the time-scales of the domain formation pro- 
Artificial fl] 12 [3) as well biological membranes fl] g] [5] C ess. Studies on the kinetics of phase transitions in artificial 
display melting transitions. In their vicinity fluctuations in membranes mainly use temperature and pressure jump tech- 
lipid area, volume, enthalpy and concentration are enhanced, niques |[24] |25] [26l [27] EH] EH] [30] [3TJ , but also volume pertur- 
As a result domains form, the elastic constants are increased bation (32] [33] [34] [35] ED, pH or ion concentration changes 
EH and the permeability is higher |H] M \M [03 ■ A whole [F37], ultrasonic measurements |38] |39] and ac-calorimetry 
variety of molecules like proteins, peptides or other small gQ\ were applied. The number of relaxation processes in 
molecules influence melting transitions. Nowadays biolog- different studies varies from one to five. Relaxation times 
ical membranes are seen as being rather heterogenous [12] i n single-lipid membranes were reported to lie in the range 
instead of being homogenous fluids as suggested in the fluid- f rom ns to min ED, depending on the experiment and the 
mosaic model EH . Heterogeneities of biological membranes detection method. Most studies agreed on a slowing down of 
due to lipid domain formation have been postulated for many relaxation processes in the transition regime ll24l l26l [38l l27l 
years HI] . In the biology community, however, domain for- [28] 0J] [33] [31] [39] [34] • It was also pointed out that maximum 
mation has only gained a stronger interest since the discus- relaxation times were observed at the transition midpoint. In 
sion about 'rafts' fl3][l6]. Rafts are thought to be involved [fJTJ relaxation processes were connected to macroscopical 
into trafficking processes [17]. Domain formation in general fluctuations and it was shown that the relaxation times of the 
is believed to trigger biochemical reaction cascades and to cooperative processes and the heat capacity are connected in 
influence enzyme activity ([18] [19] [20] [21] [22] . a simple manner. 

Domains can form as a consequence of melting transi- ln earlier studies, an influence of cholesterol, dibucaine 

tions. Already in the 1970s melting transitions were mea- and peptides was found ED EQ] [36] EQ, but not systemati- 

sured in artificial J2E and in biological membranes EES cally investigated. Molecules like general anesthetics, neu- 

using differential scanning calorimetry. The increased fluctu- rotransmitters and antibiotics display a functional role in bi- 
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ological cells. The antibiotic gramicidin A is a hydropho- 
bic, channel-forming peptide [42 1. Dimers of this peptide 
induce channels of an outer diameter of about 5 A that is 
mainly permeable for monovalent cations. Neurotransmitters 
are molecules which mainly occur in the nerve system, but 
can also be found in other parts of the body. They influence 
nerve pulse propagation and are either inhibitory or excita- 
tory. An example of an excitory neurotransmitter is sero- 
tonin (5-hydroxytryptamine). Neurotransmitters have been 
discussed to also act as anesthetics l43l . Anesthesia is the 
state when pain, consciousness or other sensations are blocked. 
General anesthetics such as 1-octanol or halothane (2-bromo- 
2-chloro-l,l,l-trifiuoroethane) lead to a reversible complete 
loss of consciousness and sensation. The action of general 
anesthetics is still not fully known. Theories favoring an di- 
rect influence on protein function are at the moment favored 
0441 1451 l46l . but lipid membrane mediated mechanisms are 
also discussed. Ueda and Yoshida [47 1 claim that anesthet- 
ics action results from effects on both proteins and lipids and 
they consider the lipid/water interface. Cantor |48l l49l re- 
lates anesthetic function to the influence of anesthetics on the 
lateral pressure profile of lipid membranes. Heimburg and 
Jackson ll50ll5TI attribute anesthetic action to the influence of 
anesthetics on the lipid membrane state. 

Peptides, neurotransmitters and anesthetics all display sig- 
nificant influence on lipid melting transition. This immedi- 
ately raises questions about these molecules on the macro- 
scopic membrane properties. In this paper, we focussed on 
the influence of anesthetics, neurotransmitters and antibiotics 
on the time scales of domain formation processes. This is 
an extension of a previous study by Grabitz and collaborators 
ll3~Tl . Relaxation processes are related to the cooperative fluc- 
tuations. More specifically, it was shown that heat capacity 
and relaxation times in pure lipid membranes are proportional 
functions. Here, the study is extended towards lipid systems 
with the incorporation of drugs. It will be shown that the ad- 
dition of these molecules does not change the linear relation 
between heat capacity and relaxation times, but systemati- 
cally alter the details of the relaxation process. We discuss the 
biological importance of our findings in relation to the influ- 
ence of domains on biochemical reaction cascades, enzyme 
activity, sorting and trafficking in biological membranes. 

Materials and Methods 

Sample preparation: Lipids were purchased from Avanti 
Polar Lipids (Birmingham/AL, USA) and used without fur- 
ther purification. All lipid samples were dissolved in a 10 
mM Hepes buffer with 1 mM EDTA at pH 7. Samples 
of multilamellar vesicles (MLVs) were prepared by adding 
buffer to the lipid powder and stirring the solution above the 
main phase transition temperature using a magnetic stirrer for 
at least one hour. During this time the lipid solution was 
vortexed at least three times. Suspensions of large unilamel- 



lar vesicles (LUVs) were prepared using an Avestin extruder 
system (Avestin Europe GmbH, Mannheim, Germany). Sus- 
pensions of MLVs were extruded at 320.9 K using a filter 
with a pore size of lOOnm. DMPC/ neurotransmitter sys- 
tems were made by adding serotonin (Hydrochloride; Sigma- 
Aldrich Inc., St. Louis/MO, USA) to the buffer before dis- 
solving the lipids in the buffer/neurotransmitter solution. In 
the case of DMPC/anesthetics solutions 1-octanol (Sigma- 
Aldrich Inc., St. Louis/MO, USA) or halothane (2-Bromo-2- 
chloro-l,l,l-trifluoroethane, Sigma-Aldrich Inc., St. Louis/- 
MO, USA) was added to the already prepared suspension of 
MLVs. The solution was stirred for another 30 rain. DMPC/- 
peptide membrane systems were prepared by dissolving the 
DMPC lipids and gramicidin A as powder in organic solvent 
(dichloromethane: methanol 2:1). After mixing both solu- 
tions the lipid/gramicidin A solution was dried through heat- 
ing and a light nitrogen or air stream. The sample was kept in 
a high vacuum desiccator over night. The rest of the prepara- 
tion equaled the one of preparing suspensions of MLVs. 

Calorimetry: Differential scanning calorimetry (DSC) 
measurements were performed with a VP-DSC from Micro- 
cal (Northhampton/MA, USA) using high feedback mode at 
a scan rate of 5K/hr if not indicated otherwise. A concen- 
tration of 10 mM was used to measure the excess heat ca- 
pacity profile of a sample directly in the calorimeter cell and 
determining the transition enthalpy. In other cases sample so- 
lutions were filled into self-built pressure cells (see below). 
In these cases concentrations of either 40 mM or 50 mM 
were prepared. The capillary volume was not exactly known, 
but with the knowledge of the transition enthalpy we could 
determine absolute excess heat capacity values. Heat capac- 
ity profiles of DMPC/gramicidin A mixtures were determined 
with scan rates of 1 K/hr. Due to hysteresis effects the scan 
rate has an influence on the temperature of the maximum heat 
capacity. Curves were corrected correspondingly. Sample so- 
lutions and buffer solutions were degassed for at least 15 mm 
before the calorimetric experiments to avoid measurement ar- 
tifacts arising from possibly evolving gas bubbles. 

Pressure perturbation calorimetry: The same calorime- 
ter was used for pressure perturbation calorimetry (PPC) ex- 
periments in the isothermal mode of the calorimeter. Mea- 
surements on solutions of MLVs were performed as previ- 
ously described [31]. The sample solution was filled into a 
self-built pressure capillary which can hold pressures up to 
500 bar. This capillary was inserted into the sample cell of 
the calorimeter. The melting transition regime shifts upon 
addition of hydrostatic pressure by about lK/40bar (see ll52l 
and with citations therein). This is used to perturb the sys- 
tem. Pressure is controlled through addition or release of high 
pressure nitrogen by opening and closing manual or solenoid 
pressure valves (Nova Swiss, Effretikon, Switzerland). The 
change of pressure occurs on a time scale of around 100 ms. 
Changing the pressure on the lipid dispersions by zLAObar and 
an appropriate choice of temperature allowed to jump to dif- 
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Figure 1 : Left: Representative calorimeter responses after a negative pressure jump. The curves were obtained 
after jumps to different temperatures within the lipid melting transition ofDMPC vesicles. The further one jumps 
through the transition (for negative pressure jumps the higher the temperature) the bigger is the integrated heat 
absorption of the lipid membrane (left, upper panel,). After normalization of the curves it can be seen that they 
decay on different time scales (left, lower panel,). Right: Representative temperature change after a positive 
pressure jump recorded by a thin thermocouple. The total response is the sum of the response of the water filled 
capillary (instrument response) and the response from the lipid sample (lipid response). 



ferent points of the phase coexistence regime where domains 
started to form. Jumps into the transition regime were either 
from below the transition (—40 bar; negative pressure jump) 
or from above the transition (40 bar; positive pressure jump)e. 
In the first case the signal is endothermic, while in the other 
case it is exothermic. After a pressure change the calorimeter 
needs to compensate the heat absorption or release to keep 
the sample temperature constant. This response contains in- 
formation about the relaxation times of the domain formation 
process. Representative calorimetric responses after negative 
pressure jumps (—40 bar) are shown in the left panel of fig.Q] 
(top). The curves are results after jumps to different points of 
the transition defined by different sample temperatures. The 
area below the curves, i.e. the total heat absorption of the 
lipid membranes during the equilibration process increases 
the further one jumps through the transition (here: increasing 
temperature). The normalized signals show that the relax- 
ation process occurs on different time scales (fig.[H bottom). 
A sample cell containing buffer solution also results in a char- 
acteristic response. This is the response function of the instru- 
ment. The observed calorimetric signal from the equilibration 
process is the convolution of the real signal (called the "sam- 
ple response") with the instrument response function. The 
response after a pressure jump always contained one small 
component that contain the perturbation of the pressure cell 



itself and small contributions from lipid membrane processes 
faster than the resolution of the instrument. The dominant 
part of the signal consists of a single exponentially decaying 
heat absorption or release from the lipid membranes after a 
perturbation. Faster relaxation processes only contribute to a 
minor degree to the overall signal (see discussion concerning 
Fig. 6). The above method allows a time resolution of about 
3 -4s. 

The relaxation times of lipid dispersions are related to 
the heat capacity. Broader melting profiles display smaller 
heat capacity values in the transition. Under these conditions 
faster relaxation processes are found (Grabitz2002) for which 
the resolution of 3 — 4 s may not be good enough. For this 
reason, a refined experimental setup with increased time res- 
olution was used especially for measurements on LUV dis- 
persion that display melting peaks that are about 10 times 
broader than the MLV transitions. We only performed pos- 
itive pressure jumps. The pressure release or addition was 
controlled by two computer controlled solenoid valves (Nova 
Swiss, Effretikon, Switzerland) and the time scale of the pres- 
sure relaxation was faster than 90 ms. The pressure cell con- 
sisted of two capillaries instead of only one. A K-thermocou- 
ple with a grounded hot junction of a diameter of 0.25 mm 
(Thermocoax, Stapelfeld, Germany) was put into the sam- 
ple capillary of the pressure cell. Relaxation times were de- 
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termined from a direct measurement of changes in temper- 
ature of the sample solution after the pressure jump. The 
thermocouple signal was amplified by a Nanovolt Preamp 
Model 1801 (Keithley Instruments Inc., Cleveland/OH, USA) 
and recorded by a Keithley Multimeter 2001 (Keithley Instru- 
ments Inc., Cleveland/OH, USA). 

As in the case analyzing the calorimeter response the to- 
tal temperature change is seen as a convolution of the "instru- 
ment response" and the "sample response". This is displayed 
in the right panel of fig. [T] The total temperature change 
(solid curve) after a positive pressure jump can be divided 
into a change due to the "instrument response" (dotted curve) 
and the "lipid response" (dashed curve). The latter one is 
again modeled by a convolution of a single exponentially de- 
caying heat release of the lipid membrane and the instrument 
response. For details we refer to IT3T1 l53l . The contribution 
from the sample cell (only filled with water) was subtracted 
from the total temperature change and only the remaining sig- 
nal was analyzed. 

Pressure perturbation calorimetry experiments were al- 
ways conducted on a few successive days. Melting profiles 
were controlled after each measuring day. During these days 
heat capacity curves sometimes broadened slightly or shifted 
towards higher or lower temperatures. The shift could be eas- 
ily corrected later. Below it will become obvious that a broad- 
ening of the curve means that relaxation times are influenced. 
The broadening, however, was small so that the determination 
of relaxation times was in the range of experimental error. 

Theory 

In an earlier paper from our laboratory a theory of the relax- 
ation times of lipids in the melting regime was derived on 
the grounds of non-equilibrium thermodynamics B3Tj. This 
theory made use of the experimental fact that enthalpy, vol- 
ume and area are proportional functions in this temperature 
regime. Therefore, perturbation of the lipid samples by tem- 
perature, bulk or lateral pressure results in the same relaxation 
process. It was found that one expects a proportional relation 
between relaxation time and excess heat capacity. 

The theory is based on the assumption that the distribution 
of enthalpy states at a given temperature can be described by 
a Gaussian distribution. This is correct when the system is a 
continuous transition that is neither of first order nature nor at 
a critical point: 
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where H is the enthalpy, (H) is the mean enthalpy and a 2 = 
(H 2 ) — (H) 2 is the variance. The Gibb's free energy depends 
logarithmical on this distribution |54l : G = RT In P. Using 
this one finds that the entropy S can be approximated as a 
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where R is the gas constant. 

Throughout a melting transition changes in area, volume 
and enthalpy are proportional and, therefore, one only finds 
one independent fluctuation given by a = H — (H) . The ther- 
modynamic force driving the lipid sample back to equilibrium 
(after a perturbation) is given by X = dS/d(H — (H)), and 
the flux of heat by J = d(H - (H))/dt. This is in fact the 
observable in our calorimetric experiments. According to On- 
sager [55 56] the flux is proportional to the thermodynamics 
force. Therefore 
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using the phenomenological constant L. This leads to a sim- 
ple differential equation that is solved by an exponential de- 
cay 
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introducing the relaxation time t. Using the identity c p = 
(H - (H)) 2 /RT 2 = <j 2 /RT 2 (fluctuation theorem) one ar- 
rives at 

T 2 



L Cp ' 
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where T is the temperature and Lisa phenomenological con- 
stant. The constant L has to be obtained from experiments. 
For the details of the derivation we refer to PP . Note, how- 
ever, that in the cited paper a factor of RT was omitted in the 
derivation so that an unit error occurred. This, however, does 
not change the message of the previous paper that relaxation 
times and heat capacity are proportional functions. 

Results 

The non-equilibrium thermodynamics theory briefly outlined 
above predicts a proportionality between heat capacity and re- 
laxation times for the cooperative processes with large excess 
heat capacity. Important for this prediction is the cooperative 
nature of the process rather than the chemical composition of 
the membrane. In ||3T| it is shown for different one compo- 
nent lipid systems (MLV and LUV) and a DMPC/cholesterol 
mixture (MLV) that one finds a linear relationship between 
heat capacity and relaxation times in pressure perturbation 
calorimetry. Maximum relaxation times are up to about a 
minute and they are decreased by the addition of cholesterol 
in a manner closely related to the effect of cholesterol on the 
heat capacity profile. In the present study we extended the 
range of the membrane systems. In particular, we investigated 
the influence of small drugs on domain formation processes 
and the related relaxation times. For this it was necessary 



4 



250x10 



Cm 



200 



150 



100- 



50- 



50 mM DMPC + 
1-octanol: 

: 0.0 wt% 

: 0.1 wt% 
— -:0.5wt% 




50 mM DMPC + 
serotonin: / 



0.0 mM 
0.5 mM 
3.0 mM 




50 mM DMPC + 
gramicidin A: 

:0.0mol% 

: 0.5 mol% 

: 1.5 mol% 



M.T.f ■*■■•■.. . 




296.0 



296.5 



297.0 296.0 



296.5 



297.0 296.0 



296.5 



297.0 



temperature [K] temperature [K] temperature [K] 

Figure 2: Heat capacity profiles of various multilamellar lipid dispersions in the presence of drugs. Left: DMPC 
plus various concentrations of 1-octanol. Center: DMPC plus different concentrations of serotonin Right: DMPC 
plus different concentrations of gramicidin A. In all three panels the drugs broaden the lipid melting profile and 
shift it to lower temperatures. 



to improve the time resolution of our calorimetric relaxation 
measurements. The latter part is important for large unilamel- 
lar membrane systems that display smaller heat capacities and 
broader transitions than multilamellar preparations. 

First, we determined the influence of molecules such as 
general anesthetics, neurotransmitters and antibiotics on the 
melting behavior of lipid membranes by differential scanning 
calorimetry. In fig. |2] the heat capacity profiles of 50 mM 
multilamellar DMPC vesicles with different concentrations 
of 1-octanol (left panel), serotonin (center panel) and grami- 
cidin A (right panel) are displayed. The transition enthalpy of 
these systems was determined to be 21.5 kJ/mol. The gen- 
eral finding is that the addition of the respective molecules 
broadened the transition profile and shifted it to lower tem- 
peratures. The total melting enthalpy remains unaltered. In 
the following we wanted to understand in which way relax- 
ation processes were influenced in the respective systems. 

In fig. [3] heat capacity profiles and relaxation times of 
a 50 mM DMPC (MLV) solution are displayed. The heat 
capacity curve shown is the same as in fig. Relaxation 
times were determined from measurements using the same 
lipid dispersion as for the DSC scans. The melting transi- 
tion profile of lipid membranes with increased pressure was 
shifted to higher temperatures, whereas the shape of the pro- 
file was not influenced. Therefore, the relaxation time profile 
obtained from positive pressure jumps were temperature cor- 
rected with AT = —0.9 K (open squares). Results from 
negative pressure jumps are given by open circles. The time 



resolution of the setup is indicated by a dashed line. Relax- 
ation times faster than this experimental response time cannot 
be determined accurately. 
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Figure 3: Heat capacity curve ( solid curve) of a 50 mM 
multilamellar DMPC sample and relaxation times ( open 
circles and squares) obtained from pressure perturba- 
tion experiments on the same sample as measured by the 
calorimeter. The different markers refer to different ex- 
perimental situations where pressure is released (open 
circles) or added (open squares). The change in relax- 
ation times correlates nicely with the heat capacity pro- 
file. The dashed line indicates the lower resolution limit 
of the experiment (3-4 seconds). 
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Comparison of heat capacity and relaxation times shows 
the linear relation of these two functions. In the proximity of 
the melting transition relaxation processes slows down. The 
maximum relaxation time measured is about (31.0 ± 1.5) s. 
With a maximum heat capacity value of 232 kJ/mol K and 
the temperature at the transition midpoint of 296.77 K we 
found a phenomenological constant (see eq.[5j of (6.6±0.5) • 
10 8 J K/mol s. The previous experiment was performed 
while recording the heat compensation of the calorimeter (with 
a time resolution of about 4 seconds). Performing experi- 
ments on a new suspension of 40 mM DMPC (MLV) and 
evaluating the temperature change with the thermocouple in- 
serted into the pressure cell (with a time resolution of about 
0.3 seconds) after a positive pressure jump of + I5bar yielded 
a phenomenological constant of (6.8 ± 0.3) • 10 s J K/mol s. 
Both phenomenological constants obtained agree within er- 
ror with each other. This demonstrates that our two different 
setups record the same time scales. 

In the following we investigated in which way other small 
molecules and peptides influence the relaxation behavior of 
lipid membranes. Therefore, we performed experiments ad- 
ding different concentrations of the anesthetic 1-octanol, the 
neurotransmitter serotonin and the antibiotic gramicidin A to 
the DMPC lipid membranes. As seen in fig.[2]these molecules 
broadened the melting transition regime and shifted it to lower 
temperatures. This has also an influence on the relaxation be- 
havior of the lipid membranes. In fig. 0] measurements on the 
three different systems are displayed. Heat capacity profiles 
are taken from fig. [2] Results were obtained from analyz- 
ing the calorimetric response after jumps of ±40bar. Results 
from negative and positive pressure jumps are not indicated 
by different symbols. The dashed line always indicates the 
time resolution of the setup. Both heat capacity and relax- 
ation time scale of the different panels are different. In all 
cases we found a proportionality between heat capacity and 
relaxation times. Phenomenological constants were calcu- 
lated to be (6.3 ± 0.3) • 10 8 J K/mol s (DMPC plus 0.1 
wt % 1-octanol), (6.5 ± 0.5) • 10 8 J K/mol s (DMPC plus 
0.5 mM serotonin) and (5.8 ± 0.5) • 10 8 J K/mol s (DMPC 
plus 1.5 mM gramicidin A), i.e. within error L is independent 
of the system. Note that deviations close to 3-4 seconds are 
most likely a consequence of the finite time resolution of the 
experiment. 

A series of PPC experiments were performed using dif- 
ferent concentrations of the respective molecules (data not 
shown). The maximum relaxation times obtained and the 
phenomenological constants determined are displayed in ta- 
ble . As it is seen in fig. E] the addition of the chosen small 
drugs led to a broadening of the heat capacity curves and to a 
decrease in the maximum heat capacity value. The maximum 
relaxation time was influenced in the same way as is the max- 
imum heat capacity value. The whole relaxation time profile 
behaves the same as does the heat capacity curve. The addi- 
tion of molecules alter the relaxation behavior of lipid domain 
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Figure 4: Comparisons between heat capacity profiles 
(solid curves) and relaxation times (open circles) of dif- 
ferent multilamellar DMPC systems. Top: 50 mM 
DMPC plus 0.1 % 1-octanol. Center: 50 mM DMPC 
plus 0.5 mM serotonin. Bottom: 50 mM DMPC plus 
1.5 mM gramicidin. The drugs alter the temperature 
dependence of the heat capacity and the magnitude and 
temperature dependence of the relaxation times. The 
dashed line indicates the resolution limit of the experi- 
ment (3-4 seconds). 



formation in a systematic way. 

Further experiments were conducted on solutions with lar- 
ge unilamellar vesicles. We investigated LUVs of DMPC 
lipids with and without addition of the anesthetics halothane. 
We also studied LUVs of DPPC. The respective heat capacity 
values and relaxation times are plotted in fig. [5] Relaxation 
times were obtained using a thermocouple inserted into the 
pressure cell using positive pressure jumps. The dashed line 
represents the time resolution of the setup. In all cases we 
found again that heat capacity and relaxation times are pro- 
portional to each other. Halothane has a similar effect on the 
melting behavior as the anesthetic octanol. It broadened and 
shifted the heat capacity curve to lower temperatures. The 
maximum relaxation time was again decreased compared to 
that of the DMPC LUVs without halothane. The phenomeno- 
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Figure 5: Heat capacity profiles (solid lines) and relax- 
ation times (symbols) of three different large unilamellar 
vesicle systems. Relaxation times were recorded with a 
thin thermocouple. Top 40 mM DMPC. Center: 40 mM 
DMPC plus 20 wt % halothane. Bottom: 40 mM DPPC. 
Again a proportionality between heat capacity and re- 
laxation times was found. Halothane acts similar on the 
lipid membrane as 1-octanol. The dashed line indicates 
the resolution limit of the experiment (around 350 mil- 
liseconds). 



to test whether there are contributions from faster relaxation 
processes we analyzed the heat absorption after negative pres- 
sure jumps as a function of temperature in comparison to the 
transition enthalpy obtained from scanning calorimetry. The 
transition enthalpy of the DMPC (MLV) solution (see the heat 
capacity profile in fig.© is shown as a solid curve in the left 
and right panel of fig.|6] The total heat absorption is displayed 
in the left panel of fig. [6] It is proportional to the transition 
enthalpy. All relaxation data contained contributions of a fast 
component with a time scale corresponding to the instrument 
response time. The heat absorption due to this component 
also showed a proportionality to the transition enthalpy (right 
panel of fig. |6). This means that at least one faster relax- 
ation process is present which occurs on a time scale faster 
than our experimental time resolution. Analyzing all mea- 
surements we found that the contribution of faster relaxation 
processes is lower than 10 % of the total heat absorption. Us- 
ing the thermocouple signal as a detection method we did not 
detect faster relaxation processes. 



Lipid Membrane System 


Tm/s 




DMPC (MLV) 


31.0 ± 1.5 


6.6 ±0.5 


+ 0.1 wt% 1-octanol (MLV) 


17.7 ±0.4 


6.3 ±0.3 


+ 0.5 wt% 1-octanol (MLV) 


14.8 ±0.2 


6.3 ±0.3 


+ 0.5 mM serotonin (MLV) 


21.6 ± 1.2 


6.5 ±0.5 


+ 3.0 mM serotonin (MLV) 


9.0 ±0.9 


7.5 ±0.8 


+ 0.5 mol% gramicidin A (MLV) 


20.8 ±0.8 


7.9 ±0.5 


+ 1.5 mol% gramicidin A (MLV) 


13.1 ±0.8 


5.8 ±0.5 


DMPC (MLV) 


25.2 ±0.1 


6.8 ±0.3 


DMPC (LUV) 


1.5 ±0.2 


15.7 ±3.0 


+ 20 wt% halothane (LUV) 


1.0 ±0.1 


20.0 ± 1.0 


DPPC (LUV) 


2.0 ±0.1 


18.4 ± 1.0 



Table 1 : Given are the relaxation times and phenomeno- 
logical constants at the transition midpoint of the studied 
systems. In the upper parts relaxation times were deter- 
mined from the calorimeter response. In the lower part 
they were obtained from the temperature change curve in 
the sample capillary. 



logical constants for the LUV systems are somewhat higher 
than those of the MLVs (see table |5). The phenomenolog- 
ical constants obtained from measurements on MLVs sug- 
gest that there is a constant independent of the lipid mem- 
brane composition. The same is true for LUVs, while we cal- 
culated different phenomenological constants for MLVs and 
LUVs. From our results we calculated averaged phenomeno- 
logical constants of (6.7 ± 0.7) • 10 8 J K/mol s (MLVs) and 
(18.0 ± 2.2) • 10 s J K/mol s (LUVs). The reason for this is 
unknown. One can speculate on whether relaxation processes 
are influenced by volume changes of the vesicles in the case 
of LUV or by bilayer-bilayer interactions in the case of MLV. 

In the literature up to five relaxation processes in the time 
range from ns to s were reported for one lipid component 
systems ATI . In this study we assumed only one relaxation 
time, which describes domain formation processes. In order 



In total we can state that domain formation processes are 
related to fluctuations in enthalpy at large excess heat capac- 
ities. Relaxation times and heat capacity are linear functions 
(within error). This linear relation is still true in the presence 
of drugs. However, the temperature dependence of both heat 
capacity and relaxation times are altered by drugs in a simple 
systematic manner. Fast relaxation processes contribute less 
than 10% to the overall heat. Thus, they probably do not rep- 
resent the fluctuations in domain size linked to the lipid chain 
state. 

Discussion 

The melting of lipid membranes is a cooperative event in- 
volving many lipid molecules. Within the melting transition 
one finds domain formation, but also changes in the elastic 
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temperature [K] temperature [K] 

Figure 6: The solid curve in the left and right panel show the transition enthalpy ofDMPC multilamellar vesicles 
as obtained from integrating the excess heat capacity profile. The total heat compensation of the calorimeter after 
a pressure jump from below the transition to a given point within the transition is proportional to the transition 
enthalpy (circles, left panel). The signal typically contains one single-exponential slow component and a minor 
contribution from fast processes that show up as a signal with the characteristic time constant of the calorimetric 
setup. The contribution from the fast processes ( circles, right panel ) also displays a proportionality to the integrated 
heat. This means that there is at least one relaxation processes faster than the time resolution of our experiment. 
Its contribution, however, is less than 8% of the overall signal. Thus, the relaxation process of the cooperative 
events is well approximated by a single relaxation process. 



constants. If pressure, temperature or any other relevant in- 
tensive thermodynamic variable is slightly changed, the lipid 
system changes its state and domains grow or shrink while 
absorbing or releasing significant amounts of heat. Thus, the 
domain formation process can be followed by monitoring the 
heat change. If we talk about relaxation times within the lipid 
melting transitions we therefore basically refer to the time 
scale of domain growth. These processes have been studied 
by some groups before using calorimetry ll33l[34ll35l[36ll3Tl . 
One typically finds slow relaxation in transitions. This phe- 
nomenon is known as 'critical slowing-down'. 

Studies of relaxation processes in artificial membranes by 
use of pressure perturbation calorimetry have been performed 
by our group before ||3T1 . We used the same and a refined ver- 
sion of the setup as described in ||3T1 to extend the previously 
published studies. Relaxation processes in LUVs and the in- 
fluence of small drugs on the relaxation behavior of lipid do- 
main formation in LUVs and MLVs were investigated. The 
typical time scales of domain formation processes are related 
to cooperative fluctuations in enthalpy and in all systems we 
found a proportionality between heat capacity and relaxation 
times. The relaxation processes depend on the final state 
rather than on the initial state. A remarkable finding is that the 
drugs used modulate the time-scale of domain formation pro- 
cesses in a similar manner as they influence the heat capacity 



profiles. They tend to broaden both heat capacity profiles and 
relaxation time profiles. After addition of anesthetics, neuro- 
transmitters and antibiotics the maximum relaxation time is 
decreased in the same way as is the maximum heat capacity 
value . Therefore, these drugs alter static and kinetic aspects 
of domain formation. 

Our findings support the theoretical prediction from Gra- 
bitz et al. |[3l1 that heat capacity and relaxation time display 
a linear relation not only for pure lipid membranes, but also 
for membrane mixtures and mixtures with drugs that asso- 
ciate with membranes. The proportional constant between re- 
laxation time and heat capacity takes the form T 2 /L, where 
L is an phenomenological coefficient originating from On- 
sager's phenomenological equations. We calculated average 
phenomenological constants to be (6. 7 ±0.7) • 10 8 J K/mols 
for MLVs and (18.0 ± 2.2) • 10 s JK /mol s for LUVs made of 
DMPC and DPPC. The constant L is different for MLV and 
LUV by about a factor of 2.5. Surprisingly, however, the addi- 
tion of drugs does not change the phenomenological constant 
for a given vesicular system even though they have signifi- 
cant influence on the heat capacity profiles. With measure- 
ments on suspensions of multilamellar DMPC vesicles we 
obtained phenomenological constants of 6.6 • 10 s J Kj mol s 
and 6.8 • 10 s J K/mol s using the two different detection 
methods. They agree within error. Thus, the difference in re- 
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suit does not seem to be related to the method of detection. 
We do not know the origin the difference between MLV and 
LUV. However, MLV display much more cooperative transi- 
tions because of coupling between adjacent membranes. The 
fluctuations therefore contain a coupling in the third dimen- 
sion that might cause small differences in the relaxation be- 
havior. The difference may also be related to the change of 
vesicle size while changing the lipid membrane state and area. 
This may involve permeation of water into the vesicles. 

Anesthetics, neurotransmitters and antibiotics influence 
the melting behavior of lipid membranes. The anesthetic oc- 
tanol shifts profiles to lower temperatures. Long chain alco- 
hols (with chain length larger than 10) shift them to higher 
temperatures lf57l . The latter ones do not display anesthetic 
potency. Other anesthetics like halothane (this study and un- 
published data) or methoxyfiurane 1 58 ] also lead to a decrease 
of the melting temperature. It has been stated that clinically 
relevant concentrations show effects on the shift of melting 
temperatures ll58l . The effective potency of an anesthetic 
is correlated to its ability to deplete melting transition tem- 
peratures ll59l . Serotonin possibly also acts as an anesthetic 
||431 and influences the phase melting transition. This has 
also been shown for other neurotransmitters, like GABA (7- 
aminobutyric acid), dopamine, caffeine and fluoxetin [60|. 
The interactions with peptides or proteins and other small 
molecules like anesthetics and neurotransmitters influence the 
domain structuring of lipid membranes in the phase melting 
transition regime ll6D . All these molecules influence relax- 
ation processes in a systematical way. 

In previous studies kinetic aspects of melting transitions 
in artificial single-component lipid membranes have been stud- 
ied with a whole variety of different techniques. This study, 
however, is the first one systematically investigating the ef- 
fects of various small drugs on relaxation processes in single 
lipid membranes and relating it to the enhancement of fluctua- 
tions in the transition regime. Van Osdol and co-workers 11361 
studied the effect of the anesthetics dibucaine on relaxation 
processes. They claimed that even though dibucaine clearly 
shows an influence on the melting behavior only a small al- 
teration of relaxation processes has been found. This finding 
does not agree with ours, where we have found that the drugs 
influence relaxation processes similar to their influence on the 
heat capacity. 

In other studies relaxation times faster than the ones we 
measured and up to five relaxation processes were reported 
fl24l ESI |26l EH E2 [62l SI] [38] [391 . Different spectroscopic 
detection methods were used. Previously, it has been dis- 
cussed that some of the experiments might have problems 
with a well-defined temperature, due to the necessity of hav- 
ing windows for the detection. These might act as heat sinks 
|[3"T1 . In favor of longer relaxation times are estimations us- 
ing ac-calorimetry ll40l . The authors found relaxation times 
faster than 120s (MLV of DMPC) and 260s (MLV of DPPC). 
They, however, did not give relaxation times in dependence of 



temperature. As mentioned, relaxation processes have been 
reported from the ns to minute regime. One reason might lie 
in an inaccurate temperature control in many optical meth- 
ods. However, another reason, which seems likely to us is that 
the different detection methods probe different relaxation pro- 
cesses. In PiTI the different relaxation times where assigned 
to relaxation processes reflecting chain properties like the for- 
mation of kinks, free rotations of headgroups or the formation 
of complex rotational isomers and to macroscopical proper- 
ties, i.e. the formation of clusters. Our study does not exclude 
the existence of several faster relaxation processes, but the fo- 
cus was only on one relaxation process, namely the formation 
of domains which is responsible for most of the heat absorp- 
tion. We found evidence that there is at least one faster relax- 
ation process which has a minor contribution to the total heat 
absorption (less than 10%). Most previous studies agreed on 
a slowing down of relaxation processes in the melting transi- 
tion regime l24ll2^l3^I^I^CTl3lll3ll|39ll34l. 

Biological membranes also display melting transitions [4 
l23l . From the knowledge of the phenomenological constant 
and the excess heat capacity profile were estimated relaxation 
times in the melting transition regime of biological mem- 
branes [31]. The authors used a DSC measurement of bovine 
lung surfactant from Ebel et al. Q and a phenomenologi- 
cal constant similar to the one for DMPC and DPPC MLV 
vesicles. They estimated that relaxation times in biological 
membranes would lie in a time regime of up to 120 ms in 
the transition regime. We have calculated two different phe- 
nomenological constants for MLVs and LUVs but the order 
of magnitude is unaffected. The time regime of these domain 
formation processes is in a regime which seems to be impor- 
tant in biology. These relaxation times do not only give the 
typical time scales of domain formation processes, but also 
the life time of domains. Domains are subject to fluctuations 
ll6D and the smaller they are the shorter their life time. Very 
broad melting profiles typically result in small domains. In 
the ongoing discussion on 'rafts' this may be of some rele- 
vance. Their size is thought to be smaller than 100 nm lfl6ll 
and therefore should they display short life times. Additional 
to their small size the short life time may make it difficult to 
really detect them. 

Domain formation processes in general have been argued 
to be important in the control of biochemical reaction cas- 
cades in biological membranes lfl8l [191 l20l |2~T1 [22] ■ Since 
drugs like anesthetics, neurotransmitters and antibiotics also 
alter the lateral membrane structure and relaxation times they 
might show an indirect action on enzymatic processes. For 
several enzymes it has been shown that the existence of do- 
mains is necessary for their function [63, 64 65 66, 67ll68l 
[69]. In the case of phospholipase A2 is was stated that the 
magnitude of fluctuations determines the activity of the en- 
zyme [70|. It may be easier for the phospholipase to reach 
the target site at the lipid when the fluctuations are large. The 
kinetics of lateral membrane structure might generally act as 



9 



a control of biochemical reactions and their time scales. 

The action of general anesthetics is still not fully under- 
stood. Many groups presently favor a picture where anes- 
thetics directly influence the function of proteins [44, 45 , 46 1. 
Theories about the action of general anesthetics mediated by 
the lipid membrane are put forward by several other groups 
ll48l |49l l50l ISTl . The later models are in our opinion more 
convincing because the action of anesthetics is known to be 
linearly related to their solubility in lipid membranes, includ- 
ing the noble gas Xenon that is inert and therefore unlikely 
to have specific interactions with macromolecules. Think- 
ing about lipid-mediated anesthesia, a possible implication 
of our findings might be related to ion permeation through 
lipid bilayers. The predicted relaxation time of biological 
membranes in the several 10 ms range is just the typical time 
scale of the opening and closing of ion channel proteins. Re- 
cently, several groups started to argue that the lipid environ- 
ment has a direct influence on the channel opening statistics 
ll7Tll72ll73l . For this reason anesthetics and neurotransmitters 
may display an indirect influence on these proteins via their 
action on the lipid membrane. This implies that the chan- 
nel lifetimes may be coupled to the relaxation time scales of 
the lipid membranes, in particular since they are obviously of 
right order. It was also reported that permeation through pure 
lipid bilayers in the absence of proteins is enhanced in the 
melting transition regime and reaches a maximum at a tem- 
perature of maximum heat capacity. This finding was related 
to the magnitude of fluctuations [IS] |9j [TDl . Measurements 
on black lipid membranes in the melting transition regime 
yielded characteristic fluctuations of the transmembrane ion 
currents ll74l l75l l76l . No such currents were found for lipid 
membranes at temperatures well above or below the tran- 
sition midpoint temperature of one component lipid mem- 
branes [75]. The life times of these lipid pores within the 
phase transition regime were found to be as long as several 
seconds, in agreement with the relaxation time scales found 
by us. 

In this study we have shown that relaxation times of do- 
main formation processes in simple model systems are pro- 
portional to the excess heat capacity. Relaxation processes 
slow down in the melting regime. Several drugs studied by us 
broaden the melting transition profiles and shift them to lower 
temperatures. This influences the relaxation behavior in a 
simple manner. This is presumably true for all drugs, pep- 
tides and proteins influencing the melting behavior of lipid 
membranes. Clearly, the free energy of membranes includes 
the chemical potentials of all membrane associated molecules 
and it is therefore just a consequence of thermodynamics that 
the functions of states change in a coherent manner due to the 
change in one of these variables. This is likely to be of great 
importance for the general understanding of the function of 
biological membranes. 
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